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ABSTRACT 

We present multi-frequency radio observations of the 2010 nova event in the symbiotic binary V407 Cygni, obtained 
with the Karl G. Jansky Very Large Array (VLA) and spanning 1-45 GHz and 17-770 days following discovery. 
This nova — the first ever detected in gamma rays — shows a radio light curve dominated by the wind of the Mira 
giant companion, rather than the nova ejecta themselves. The radio luminosity grew as the wind became increasingly 
ionized by the nova outburst, and faded as the wind was violently heated from within by the nova shock. This 
study marks the first time that this physical mechanism has been shown to dominate the radio light curve of an 
astrophysical transient. We do not observe a thermal signature from the nova ejecta or synchrotron emission from 
the shock, due to the fact that these components were hidden behind the absorbing screen of the Mira wind. We 
estimate a mass-loss rate for the Mira wind of M w ~ 1 0 6 * M 0 yr 1 . We also present the only radio detection of 
V407 Cyg before the 2010 nova, gleaned from unpublished 1993 archival VLA data, which shows that the radio 
luminosity of the Mira wind varies by a factor of >20 even in quiescence. Although V407 Cyg likely hosts a 
massive accreting white dwarf, making it a candidate progenitor system for a Type la supernova, the dense and 
radially continuous circumbinary material surrounding V407 Cyg is inconsistent with observational constraints on 
the environments of most Type la supernovae. 

Key words: binaries: symbiotic - novae, cataclysmic variables - radio continuum: stars - stars: individual 
(V407 Cygni) 
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1. INTRODUCTION 

On 2010 March 10, the symbiotic star V407 Cygni hosted 

the first nova explosion ever detected in gamma rays (Cheung 

et al. 2010; Abdo et al. 2010). This explosion also constitutes the 

second example of an “embedded” nova thoroughly observed at 

a range of wavelengths, presenting an opportunity to trace the 

interaction of a blast wave with dense circumstellar material in 

detail (the first example was RS Ophiuchi, with recent outbursts 

in 1985 and 2006; Bode 1987; Evans et al. 2008). Many parallels 
can be drawn between RS Oph and V407 Cyg (Mikolajewska 
2010; Nelson et al. 2012), and together they form a class of 
intriguing candidates for Type la supernova (SN la) progenitors 
(Munari & Renzini 1992; Hachisu 2003; Patat et al. 2011a). 

V407 Cyg is a D-type (dust-rich) symbiotic binary featuring 
a white dwarf (WD) and a Mira giant star in a wide orbit 
(Munari et al. 1990). The giant pulsates with a period of 
763 days (Meinunger 1966; Kolotilov et al. 1998, 2003), the 
longest known period for a Mira symbiotic system. There is 

some indication of a 43 year orbital period in the long-term 
optical light curve (Munari et al. 1990), corresponding to a 


binary separation of 17 AU for a WD mass of 1 M Q and a giant 
mass of 1.5 M Q . 

V407 Cyg has undergone periods of increased activity 
throughout the past several decades, observed as peaks in the 
optical light curve and brightening emission lines, particularly 
in the years 1936 and 1998 (Hoffmeister 1949; Welin 1973; 
Munari et al. 1990; Kolotilov et al. 1998, 2003; Tatarnikova 
et al. 2003a; Esipov et al. 2012). These outbursts likely mark 
an increase in the accretion rate onto the WD (Kolotilov et al. 
2003), and so do not qualify as bona fide novae, which are 
powered by a thermonuclear runaway. 

However, the 2010 outburst was a true nova explosion, 
displaying an unprecedented brightening of the optical light 
curve (reaching a //-band magnitude of ~8 mag, as opposed to 
13.1 mag in the 1998 outburst; Kolotilov et al. 2003; Nishiyama 
et al. 2010; Munari et al. 2011) and an expansion velocity of the 
ejecta of 3200 km s 1 (Buil 2010; Abdo et al. 2010). The unique 
gamma-ray signature from this explosion has been attributed to 
particles accelerated by the interaction of the nova shock with the 
dense Mira wind (Abdo et al. 2010; Lii et al. 2011; Hernanz & 
Tatischeff 2012; Sitarek & Bednarek 2012; Martin & Dubus 
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2012). The optical properties of V407 Cyg’s nova event have 
been described by Munari et al. (2011) and Shore et al. (2011, 
2012), while the X-ray evolution was studied by Shore et al. 
(2011), Orlando & Drake (2012), Schwarz et al. (2011), and 
Nelson et al. (2012). 

This paper presents complementary radio continuum obser- 
vations of the 2010 outburst, obtained between 1-45 GHz with 
the Karl G. Jansky Very Large Array (VLA; Perley et al. 2011). 
Before the 2010 outburst, no radio detections of V407 Cyg had 
ever been reported (although, see Section 7), but shortly after 
the nova explosion, V407 Cyg was strongly detected across 
the centimeter and millimeter bands (Nestoras et al. 2010; 
Gawronski et al. 2010; Pooley 2010; Bower et al. 2010; Krauss 
et al. 2010). In Section 2, we present our data and describe the 
radio light curve. In Section 3, we present a new distance es- 
timate to V407 Cyg, based on optical Nai D absorption lines. 
Section 4 considers and dismisses the possibility that the radio 
emission from V407 Cyg is dominated by expanding thermal 
ejecta expelled in the nova. In Section 5, we hypothesize that 
the radio emission is produced by the Mira wind, similar to the 
radio-emitting ionized envelopes seen around other symbiotic 
systems (Seaquist et al. 1984), but ionized by the nova outburst. 
In Section 6, we predict the synchrotron signal from V407 Cyg 
and show that it was significant but hidden behind the absorb- 
ing screen of the Mira wind. Finally, the first radio detection of 
V407 Cyg in quiescence is presented in Section 7, using archival 
VLA data from 1993. We conclude in Section 8. 

2. THE DATA 
2.1. Observational Method 

At the time of the outburst of V407 Cyg, a new correlator was 
being commissioned as part of the VLA upgrade. We therefore 
began monitoring V407 Cyg on 2010 April 3.7, approximately 
one month after discovery (2010 March 10), as some of the 
newly upgraded facility’s first observations (Krauss et al. 2010). 
During each epoch, we attempted to observed at L (1-2 GHz), 
C (4-8 GHz), K (19-26 GHz), and Q (40-50 GHz) bands, 
providing coverage across the entire VLA frequency range. 
Under program AS 1039, we obtained data roughly twice a 
week throughout 2010 April and May. The spacing of our 
epochs gradually widened, and monitoring was continued until 
2012 April via programs 10B-233, 11A-254, and 11B-170. 
Our monitoring of V407 Cyg has now bridged more than an 
entire VLA configuration cycle, beginning in its most compact 
D configuration. During the summer (June-August) of 2011, 
the VLA was in the extended A configuration, and V407 Cyg 
was spatially resolved. The details of this imaging will be 
presented in A. J. Mioduszewski et al. (2013, in preparation), 
but flux densities were difficult to derive during this time at 
some frequencies. 

In a typical epoch, we spent two hours total to observe all 
four bands, providing ~ 10-15 minutes on source at each band. 
At K and Q bands, we obtained improved pointing solutions on 
a nearby calibrator, and used J2 102+4702 for gain calibration. 
At C band, we used J2048+4310 for gain calibration, and at 
L band, J2038+5 1 19 and J2052+3635 were used. The absolute 
flux density scale and bandpasses were calibrated using 3C48. 
The earlier data (programs AS 1039 and 10B-233) were obtained 
with 256 MHz of bandwidth per band, split between two 
independently tunable basebands. Later observations (programs 
11A-254 and 1 IB- 170) provided 2 GHz of bandwidth, split 
between two 1 GHz wide basebands. 


We also utilized some early data publicly available in the 
archive under program AL733 (PI Lang), obtained on 2010 
March 27.9 and 2010 March 29.8 at L and Cbands. This program 
observed with one-hour-long blocks at each receiver band, and 
used 256 MHz of bandwidth (at L band, this bandwidth was 
split between 1.4 and 1.9 GHz, while at C band the basebands 
were set adjacent to one another at 5 GHz). 

Data reduction was carried out using standard routines in 
CASA and AIPS. Gain curves and weather tables (using the 
mean of the seasonal average and the day’s measurements) 
were applied. Basic editing and calibration were carried out 
on each epoch and receiver band separately. Typically, one 
iteration of phase-only self-calibration was applied, and the two 
basebands were imaged separately. We measured flux densities 
of V407 Cyg by fitting Gaussians to the unresolved point-like 
source, except during the A configuration (epoch 26), when 
two Gaussians were fit and summed together to represent the 
higher-frequency flux densities. Calibration errors of 5% were 
assumed at L and C bands, 10% at K band, and 20% at Q band 
( Q band was particularly plagued by poor receiver performance, 
a faint gain calibrator, and inclement weather). Measurements 
are presented in Table 1 . 

Additional radio continuum data were published by Nestoras 
et al. (2010), Gawronski et al. (2010), Pooley (2010), and 
Bower et al. (2010) in Astronomer’s Telegrams. Unfortunately, 
the facilities and frequencies used for these measurements are 
heterogeneous, so they are difficult to integrate into our VLA 
light curve and are not used in our analysis (although a spectral 
index from the Nestoras et al. data is listed in Table 2). 

2.2. Results 

Our multi-frequency radio light curve for V407 Cyg, plotted 
in Figure 1 , shows a gradual rise and decline of radio luminosity 
over the two years following the nova outburst. Due to observa- 
tional complications, there were occasionally small shifts in the 
central frequencies of our bands, as shown in Table 1 ; therefore, 
for example, the “7.4 GHz” light curve actually includes data 
centered at frequencies between 7.0 and 7.9 GHz. To minimize 
dependence on model assumptions, no spectral index correction 
was applied to our light curve plots to account for this small 
difference. 

The first VLA observations were obtained 17 days after 
outburst at only 5 GHz, and V407 Cyg is already radio bright at 
this time (~7 mJy). Overall, the light curve rises and declines, 
with the higher frequencies peaking at larger luminosities and 
declining first. At 45 GHz, V407 Cyg is already fading by the 
time of our first reliable (J-band observation on Day 33. At 
frequencies between 4.9 and 25 GHz, the source continues to 
brighten for a month or two longer, before turning over and 
declining. At 1.4 and 1.8 GHz, the flux density rises until Day 
~70, and then remains roughly constant at this peak flux density 
for another ~ 100 days before gently declining. By Day 300, the 
flux densities at all frequencies converge to ~3.5 mJy and then 
remain roughly constant for the next year, through the time of 
writing of this paper. 

The entire radio evolution of V407 Cyg is relatively gentle and 
slow, showing gradual rises and broad peaks in the 1 .4-25 GHz 
light curves. Typical expectations for novae attribute the radio 
emission to expanding thermal ejecta, with luminosity increas- 
ing in proportion to area while it is optically thick (Seaquist & 
Bode 2008; Section 4) — so, L v oc t 2 . However, the rise of the 
radio light curves in V407 Cyg is significantly slower than this, 
scaling as L v oc t°- 5 ~ L0 . 
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Table 1 

VLA Observations of V407 Cyg 


Date 

(UT) 

t — t a 
T Iq 

(days) 

Epoch 

Conf. 

Freq 

(GHz) 

S v 

(mJy) 

Freq 

(GHz) 

Sv 

(mJy) 

Freq 

(GHz) 

Sv 

(mJy) 

Freq 

(GHz) 

Sv 

(mJy) 

2010 Mar 27.5 

17.1 

i 

D 





5.0 

6.39 ± 0.35 



2010 Mar 29.8 

19.0 

2 

D 

1.5 

-0.12 ±0.72 

1.9 

1.93 ±0.53 

5.0 

7.03 ± 0.37 



2010 Apr 3.7 

23.9 

3 

D 

1.4 

2.98 ± 0.52 

1.8 

3.47 ± 0.39 

4.8 

6.88 ± 0.37 

7.0 

8.81 ±0.45 





19.0 

20.49 ± 2.05 

25.0 

24.55 ± 2.46 



49.0 

44.98 ± 22.55 

2010 Apr 8.7 

28.9 

4 

D 





4.8 

8.99 ± 0.47 

7.0 

11.88 ±0.60 





19.0 

23.65 ± 2.37 

25.0 

27.71 ± 2.77 





2010 Apr 12.7 

32.9 

5 

D 

1.4 

4.04 ± 0.36 

1.8 

4.76 ± 0.38 

4.8 

10.68 ±0.55 

7.0 

12.87 ± 0.66 





19.0 

29.13 ±2.91 

25.0 

35.28 ± 3.53 

41.0 

59.89 ± 12.01 

45.0 

65.63 ± 13.18 

2010 Apr 15.7 

35.9 

6 

D 

1.4 

3.60 ± 0.38 

1.8 

5.47 ±0.41 

4.8 

12.36 ± 0.65 

7.0 

14.67 ± 0.74 





19.0 

30.90 ± 3.09 

25.0 

37.24 ± 3.72 





2010 Apr 18.7 

38.9 

7 

D 

1.4 

4.78 ±0.38 

1.8 

5.51 ±0.45 

4.8 

12.79 ± 0.66 

7.0 

15.59 ± 0.79 





19.0 

29.15 ±2.92 

25.0 

34.35 ± 3.44 





2010 Apr 22.7 

42.9 

8 

D 

1.4 

5.48 ± 0.41 

1.8 

6.66 ± 0.47 

4.8 

13.69 ±0.71 

7.0 

15.72 ±0.79 





19.0 

28.91 ±2.90 

25.0 

30.32 ± 3.05 





2010 Apr 25.7 

45.9 

9 

D 

1.4 

5.78 ± 0.43 

1.8 

6.84 ± 0.45 

4.8 

13.96 ± 0.70 

7.0 

16.07 ± 0.81 





19.0 

30.38 ± 3.03 

25.0 

33.86 ±3.39 





2010 Apr 29.6 

49.8 

10 

D 

1.4 

5.94 ± 0.39 

1.8 

7.38 ± 0.47 

4.8 

13.67 ± 0.69 

7.0 

16.88 ± 0.85 





19.0 

33.03 ± 3.30 

25.0 

37.51 ± 3.75 





2010 May 2.6 

52.8 

11 

D 

1.4 

5.76 ± 0.40 

1.8 

7.48 ± 0.46 

4.8 

13.14 ±0.67 

7.0 

16.10 ±0.81 





19.0 

29.07 ± 2.91 

25.0 

34.75 ± 3.48 





2010 May 7.6 

57.8 

12 

D 

1.4 

6.44 ± 0.44 

1.8 

7.66 ±0.51 

4.8 

14.16 ±0.71 

7.0 

17.18 ± 0.86 





19.0 

27.77 ± 2.78 

25.0 

29.49 ± 2.95 





2010 May 14.5 

64.7 

13 

D 

1.4 

5.36 ± 0.39 

1.8 

6.65 ± 0.43 

4.8 

15.02 ± 0.76 

7.9 

19.18 ± 0.96 

2010 May 11.6 

61.8 



19.0 

28.25 ± 2.83 

25.0 

30.31 ± 3.03 

41.0 

51.34 ± 10.30 

45.0 

50.66 ± 10.16 

2010 May 16.6 

66.8 

14 

D 

1.4 

5.80 ± 0.40 

1.8 

7.57 ± 0.46 

4.8 

15.48 ±0.79 

7.9 

19.88 ± 1.00 





19.0 

28.75 ± 1.44 

25.0 

29.25 ± 2.93 

41.0 

44.90 ± 9.00 

45.0 

46.31 ±9.28 

2010 May 25.6 

75.8 

15 

D 

1.4 

5.85 ± 0.42 

1.8 

8.05 ± 0.48 

4.8 

17.40 ±0.88 

7.9 

21.71 ± 1.09 





19.0 

28.36 ±2.88 

25.0 

28.80 ±2.88 

41.0 

39.97 ± 8.01 

45.0 

38.40 ±7.71 

2010 May 29.6 

79.8 

16 

D 

1.4 

5.83 ±0.41 

1.8 

7.53 ± 0.47 

4.8 

17.47 ± 0.89 

7.9 

20.62 ± 1.04 





19.0 

26.72 ± 2.67 

25.0 

27.51 ± 2.75 

41.0 

37.05 ± 7.43 

45.0 

35.88 ± 7.20 

2010 Jun 15.5 

96.7 

17 

D 

1.4 

5.61 ±0.38 

1.8 

7.54 ± 0.47 

4.8 

17.26 ±0.87 

7.9 

20.19 ± 1.01 

2010 Jun 20.6 

101.8 

18 

D 

1.4 

5.76 ± 0.44 

1.8 

7.81 ±0.49 

4.8 

17.78 ± 0.90 

7.9 

19.02 ± 0.95 

2010 Jun 21.5 

102.7 

19 

D 

1.4 

5.91 ±0.41 

1.8 

7.38 ± 0.45 

4.8 

17.45 ± 0.88 

7.9 

18.20 ±0.93 





19.0 

23.04 ±2.31 

25.0 

22.47 ± 2.25 

41.0 

25.94 ± 5.23 

45.0 

28.13 ± 5.68 

2010 Jul 27.4 

138.6 

20 

D 

1.4 

6.92 ± 0.44 

1.8 

8.23 ±0.51 

4.8 

12.93 ± 0.65 

7.9 

14.57 ± 0.73 





19.0 

14.87 ± 1.49 

25.0 

13.63 ± 1.37 

41.0 

18.67 ± 3.83 

45.0 

18.61 ± 3.85 

2010 Aug 11.4 

153.6 

21 

D 

1.4 

7.92 ± 0.52 

1.8 

8.19 ±0.46 

4.8 

11.81 ±0.60 

7.9 

12.92 ± 0.65 





19.0 

13.26 ± 1.33 

25.0 

13.30 ± 1.33 

41.0 

13.39 ±2.52 

45.0 

13.51 ± 2.75 

2010 Aug 26.3 

168.5 

22 

D 

1.4 

5.49 ± 0.38 

1.8 

6.39 ± 0.40 

4.8 

10.79 ±0.55 

7.9 

11.31 ± 0.57 





19.0 

11.87 ± 1.19 

25.0 

11.20 ± 1.12 

41.0 

14.58 ±2.97 

45.0 

15.12 ±3.08 

2010 Sep 10.3 

183.5 

23 

D 

1.4 

5.61 ±0.37 

1.8 

6.34 ± 0.36 

4.8 

9.98 ±0.51 

7.9 

10.16 ±0.51 





19.0 

10.01 ± 1.00 

25.0 

9.51 ±0.95 

41.0 

9.66 ±2.11 

45.0 

9.09 ± 1.98 

2010 Sep 30.2 

203.4 

24 

DnC 

1.4 

5.58 ±0.31 

1.8 

5.53 ±0.31 

4.8 

9.17 ±0.46 

7.9 

9.10 ±0.46 





19.0 

9.86 ± 0.99 

25.0 

10.08 ± 1.01 

41.0 

14.19 ±2.91 

45.0 

12.75 ± 2.64 

2011 Mar 5.2 

359.6 

25 

B 

1.3 

4.22 ± 0.27 

1.8 

4.64 ± 0.25 

5.0 

5.00 ± 0.25 

7.4 

4.81 ± 0.24 





20.1 

4.65 ± 0.47 

25.6 

4.74 ± 0.48 

41.0 

5.30 ± 1.08 

45.0 

5.45 ± 1.11 

2011 Jun 23.4 

469.6 

26 

A 

1.3 

3.60 ± 0.20 

1.8 

4.07 ±0.21 

5.0 

3.27 ±0.17 

7.4 

2.75 ± 0.14 





20.1 

2.54 ± 0.26 

25.6 

2.62 ± 0.27 

41.0 

4.25 ± 0.87 

45.0 

4.45 ± 0.92 

2011 Nov 4.1 

603.3 

27 

D 

1.3 

3.64 ± 0.56 

1.8 

3.86 ± 0.35 

5.0 

3.54 ±0.18 

7.4 

3.49 ±0.18 





20.1 

3.67 ± 0.37 

25.6 

3.83 ± 0.39 

41.0 

4.32 ± 0.87 

45.0 

4.47 ±0.91 

2012 Jan 08.8 

669.0 

28 

DnC 

1.3 

4.69 ± 1.00 

1.8 

3.32 ± 0.55 

5.0 

3.21 ±0.17 

7.4 

3.15 ± 0.16 





20.1 

3.14 ±0.32 

25.6 

2.99 ± 0.30 

41.0 

3.96 ± 0.83 

45.0 

5.23 ± 1.11 

2012 Apr 18.4 

769.6 

29 

C 

1.3 

2.94 ± 0.33 

1.8 

3.35 ± 0.29 

5.0 

2.91 ±0.15 

7.4 

2.80 ±0.14 





20.1 

2.91 ±0.29 

25.6 

2.97 ± 0.30 

41.0 

2.88 ± 0.60 

45.0 

2.66 ± 0.59 


Note. a Time since optical discovery of the nova event (2010 March 10.8; Munari et al. 201 1). 


Even at early times, the spectral index a (defined as S v oc v“) 
of V407 Cyg is significantly shallower than the a = 2 predicted 
for optically thick thermal emission (e.g., Seaquist & Bode 
2008). The report of Nestoras et al. (2010) and the early VLA 
data from program AL733 indicate that the spectrum may be 
steepest at early times, with a > 1.0 (see Table 2). By the time 


our full VLA monitoring effort began on Day 24, the spectral 
index was measured to be a = 0.7-0. 8. 

On Day 24, the spectrum could be fit with a single power law 
of a = 0.75 between 1.4 and 49 GHz (Figure 2). This spectral 
characterization remained roughly constant through Day ~70, 
with some evidence for gentle flattening with time to a 0.6. 
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10 20 50 100 200 500 1000 

Time Since Discovery (Days) 


Figure 1. Multi-frequency radio light curve for V407 Cyg. A variety of colors denote different frequencies spanning 1 .4 — 45 GHz, as marked on the left side of the 
plot. The time of optical discovery of the nova event is 2010 March 10.8. 

(A color version of this figure is available in the online journal. 1 


Table 2 

Spectral Indices for V407 Cyg 


Epoch 

t — t a 
T Tq 

1-45 GHz 


1-9 GHz 


19^15 GHz 


O' 

Xv ( dof ) 

Of 

Xv (dof) 

a 

Xv ( dof ) 

-i b 

11.7 

1.03 ± 0.18 

2.51(3) 



0.98 ± 0.26 

3.45(2) 

2 

19.0 

1.34 ±0.29 

0 

1.34 ±0.29 

0 



3 

23.9 

0.75 ± 0.02 

0.28(5) 

0.68 ± 0.01 

0.01(2) 

0.73 ±0.10 

0.06(1) 

4 

28.9 

0.69 ± 0.01 

0.06(2) 

0.74 ±0.19 

0 

0.58 ± 0.52 

0 

5 

32.9 

0.76 ± 0.02 

0.50(6) 

0.74 ± 0.05 

0.89(2) 

0.94 ± 0.08 

0.13(2) 

6 

35.9 

0.77 ± 0.05 

2.06(4) 

0.82 ± 0.09 

3.27(2) 

0.68 ± 0.52 

0 

7 

38.9 

0.70 ± 0.03 

1.17(4) 

0.76 ± 0.05 

0.91(2) 

0.60 ± 0.52 

0 

8 

42.9 

0.61 ± 0.05 

1.42(4) 

0.66 ± 0.05 

1.29(2) 

0.17 ±0.52 

0 

9 

45.9 

0.62 ± 0.03 

0.79(4) 

0.65 ± 0.05 

1.14(2) 

0.39 ±0.51 

0 

10 

49.8 

0.64 ± 0.01 

0.17(4) 

0.64 ± 0.02 

0.27(2) 

0.46 ± 0.52 

0 

11 

52.8 

0.61 ± 0.02 

0.40(4) 

0.61 ± 0.04 

0.75(2) 

0.65 ± 0.52 

0 

12 

57.8 

0.55 ± 0.03 

1.05(4) 

0.61 ± 0.02 

0.16(2) 

0.22 ± 0.52 

0 

13 

63 

0.63 ± 0.04 

3.01(6) 

0.74 ± 0.05 

1.75(2) 

0.72 ±0.16 

0.54(2) 

14 

66.8 

0.58 ± 0.04 

3.01(6) 

0.70 ± 0.05 

1.47(2) 

0.56 ±0.17 

0.58(2) 

15 

75.8 

0.56 ± 0.06 

6.56(6) 

0.74 ± 0.07 

3.24(2) 

0.39 ±0.12 

0.30(2) 

16 

79.8 

0.55 ± 0.06 

7.05(6) 

0.73 ± 0.08 

4.33(2) 

0.37 ±0.10 

0.20(2) 

17 

96.7 

0.74 ± 0.09 

5.38(2) 

0.74 ± 0.09 

5.38(2) 



18 

101.8 

0.68 ± 0.12 

8.30(2) 

0.68 ±0.12 

8.30(2) 



19 

102.7 

0.47 ± 0.07 

9.26(6) 

0.67 ±0.11 

8.31(2) 

0.20 ±0.10 

0.21(2) 

20 

138.6 

0.27 ± 0.05 

4.73(6) 

0.43 ± 0.04 

1.17(2) 

0.27 ±0.21 

0.80(2) 

21 

153.6 

0.19 ±0.04 

2.38(6) 

0.30 ± 0.03 

0.54(2) 

-0.03 ± 0.05 

0.05(2) 

22 

168.5 

0.28 ± 0.05 

4.88(6) 

0.43 ± 0.07 

2.89(2) 

0.27 ±0.17 

0.54(2) 

23 

183.5 

0.20 ± 0.05 

5.41(6) 

0.35 ± 0.06 

2.62(2) 

-0.09 ± 0.05 

0.04(2) 

24 

203.4 

0.25 ± 0.04 

2.82(6) 

0.33 ± 0.07 

3.85(2) 

0.36 ±0.12 

0.31(2) 

25 

359.6 

0.04 ± 0.02 

0.56(6) 

0.07 ± 0.03 

0.82(2) 

0.19 ±0.03 

0.02(2) 

26 

469.6 

-0.13 ± 0.05 

3.13(6) 

-0.17 ±0.07 

3.86(2) 

0.71 ±0.19 

0.59(2) 

27 

603.3 

0.02 ± 0.03 

0.45(6) 

-0.05 ± 0.02 

0.10(2) 

0.24 ± 0.02 

0.01(2) 

28 

669.0 

-0.01 ± 0.05 

1.15(6) 

-0.10 ± 0.07 

0.62(2) 

0.47 ± 0.26 

1.04(2) 

29 

769.6 

-0.03 ± 0.02 

0.40(6) 

-0.07 ± 0.05 

0.61(2) 

-0.07 ± 0.24 

0.06(2) 


Notes. 

a Time since optical F-band maximum (2010 March 10.8; Munari et al. 201 1). 
b From 2010 March 22/23, as documented in Nestoras et al. (2010). 
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Figure 2. Selected radio spectra showing the evolution of V407 Cyg. The 1^-5 GHz power-law fits described in Table 2 are marked as dashed red lines. In the 
bottom-right panel, comparison spectra are shown for three physically significant spectral indices. 

(A color version of this figure is available in the online journal.) 

After Day 70, the 25 GHz light curve clearly turns over and 
the spectra show flattening above ~20 GHz. By Day 100, the 
spectral flattening has cascaded down to ~6 GHz, and by Day 
200 the entire spectrum is consistent with a flatter power law of 
a = 0.25. At late times (> 350 days after outburst), the spectral 
index is even flatter, consistent with an optically thin thermal 
spectrum (a = —0.1; Figure 2). 

3. DISTANCE TO V407 CYG 

Our interpretation of the radio light curve depends on the 
distance to V407 Cyg, which is poorly known. Estimates in 
the literature range from 1.7 kpc (Kolotilov et al. 2003) to 
2.7 kpc (Munari et al. 1990), all using the period-luminosity 
relation for Miras. However, the period-luminosity relation is 
not constrained for the very long period of V407 Cyg (e.g., 

Whitelock et al. 2008), and the relation must be widely ex- 
trapolated to derive a distance. In addition, Tatarnikova et al. 

(2003b) identified lithium in the atmosphere of V407 Cyg. 

Li-rich asymptotic giant branch stars have been noted to 
lie above the period-luminosity relation, displaying anoma- 
lously high luminosities for a given period (Whitelock et al. 

2003). Therefore, distances to V407 Cyg derived from the 
period-luminosity relationship are highly uncertain and likely 
underestimate the true distance. 

To provide additional constraints on the distance using inter- 
stellar Nai D absorption lines, we revisited the echelle spectra 
published by Munari et al. (2011). Two velocity components 
were observed that were constant in time as the 2010 outburst 
evolved (Figure 3), indicating an interstellar foreground origin. 

They have heliocentric velocities of —45.7 ± l.lkms -1 and 
— 10.9 ± l.Okms -1 . In the local standard of rest (LSR), the 
corresponding velocities are —29.9 km s' 1 and +4.9 km s -1 (for 
Galactic coordinates / = 86.98, b = —0.48). 



heliocentric wavelength (A) 

Figure 3. Detail of echelle spectra around the Nai D absoiption lines. We 
plot eight distinct spectra, observed at a range of times after the nova outburst, 
with the absorption lines superimposed on the time-variable red wing of the 
He 1^.5976 line. The Nai D absorption in all spectra shows a constant two- 
component structure. 

We use the Galactic rotation curve of Brand & Blitz (1993) 
to translate the LSR velocities of the absorbing components to 
distances of 3.0 kpc and 1.5 kpc, respectively. Both components 
are strong and nearly saturated (Figure 3), indicating a signifi- 
cant amount of absorbing material. If the Milky Way’s neutral 
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interstellar medium is approximated as a slab extending 150 pc 
from the midplane, our line of sight will run through it for 
~3.5 kpc, consistent with the existence of absorption from a 
distance of 3.0 kpc. This line of sight crosses the Orion Spur 
spiral arm at a distance of ~ 1 .5 kpc and meets the large Perseus 
spiral arm at a distance of 3 kpc (Churchwell et al. 2009). 
Therefore, the interstellar Nai D absorption implies a distance 
>3.0 kpc for V407 Cyg. 

Several other authors have hypothesized that the more 
blueshifted absorption component (ulsr = —30 kins ') is due 
to the circumstellar envelope of V407 Cyg itself (Tatarnikova 
et al. 2003a; Shore et al. 2011). If this hypothesis were ac- 
curate, we would only be able to constrain the distance of 
V407 Cyg to >1.5 kpc. Tatarnikova et al. point out that the 
more blueshifted Na i D component has a very similar radial ve- 
locity as the Ca i and Li i absorption features, which presumably 
arise in the Mira’s atmosphere. However, we observe that the 
more blueshifted absorption component did not vary in strength 
relative to the less blueshifted component over the course of 
the nova outburst, implying that both absorption components 
originate in the interstellar — not the circumstellar — medium. In 
addition, while the narrow absorption component of Ha varies 
in its central wavelength throughout the nova outburst (it shifts 
blueward with time), the Nai absorption is remarkably constant 
in its central wavelength (U. Munari et al. 2013, in preparation). 
We therefore conclude that the correspondence in velocity be- 
tween the i>lsr = — 30kms _1 Nai D absorption component 
and stellar features is a coincidence, and both Na i D absorption 
components are interstellar. 

Hence, we find that V407 Cyg is located >3.0 kpc away. We 
will take a distance at the lower end of this range, D — 3.0 kpc, 
throughout this paper, as it is consistent with most other works 
on the 2010 nova, which assume a distance of 2.7 kpc. However, 
we maintain that the distance to V407 Cyg is still uncertain. 
Where appropriate throughout our work, we note how a change 
in the distance would affect our results. 

4. THERMAL NOVA EJECTA AS A POTENTIAL SOURCE 
OF RADIO EMISSION FROM V407 CYG 

Here, we compare the multi-frequency radio light curve of 
V407 Cyg with a model that has historically explained the radio 
emission from most classical novae: expanding, thermal ejecta 
spherically expelled by the nova. We show that if this mechanism 
is the dominant source of radio luminosity, the nova ejection 
in V407 Cyg must have been massive (10~ 5 to I 0“ 4 M 0 ) — in 
conflict with significant observational evidence that the ejecta 
were low in mass (10~ 7 to 10~ 6 Mq). 

In the past, most measurements of radio emission from novae 
were consistent with thermal bremsstrahlung from a ~ 1 0 4 K 
ionized sphere that is expanding in size, dropping in density, 
and gradually transitioning from optically thick to optically 
thin (e.g., Seaquist & Palimaka 1977; Hjellming et al. 1979). 
The radio emission has been interpreted as the expanding nova 
ejecta, and from basic assumptions about its density profile 
and from simple modeling of its expansion, the ejecta mass 
and energy can be estimated (Seaquist & Bode 2008). This 
scenario has mostly been applied to classical novae with minimal 
circumbinary material surrounding them, and indeed, the radio 
emission from the only other well-studied nova in a symbiotic 
system, RS Oph, has proven much more complex (Hjellming 
et al. 1986; O’Brien et al. 2006; Kantharia et al. 2007; Rupen 
et al. 2008; Sokoloski et al. 2008; Eyres et al. 2009). Still, here 
we begin with the simplest explanation for radio emission in 
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V407 Cyg, and compare our VLA light curve with expectations 
for the expanding nova ejecta alone. 

To first order, most radio light curves of classical novae can be 
fit with a simple model wherein the ejecta are instantaneously 
ejected and homologously expanding. The model requires a 
maximum and minimum velocity (u max and i.' lnm ) as input for this 
“Hubble flow,” and the ejecta are assumed to be distributed in a 
thick expanding shell with a p cx r~ 2 radial density profile. In 
this simple model, the ejecta are also assumed to be isothermal, 
with an electron temperature ~10 4 K. Additional complications 
can certainly be added to this model, like clumping, temperature 
increases, and non-spherical geometry (e.g., Hjellming 1996; 
Hey wood 2004; Roy et al. 2012). However, in most cases, the 
simplest model has been sufficient to approximate observed light 
curves. 

In Figure 4, we superimpose such models on our VLA light 
curve for V407 Cyg to test if the ejecta might account for 
the majority of the radio emission. We assume a distance of 
3.0 kpc, ii m ax = 3, 200kms _1 (Munari et al. 2011; Shore et al. 
2011), and a thick shell with u m i n /u ma x = 0.1. Although no 
Hubble-flow model provides a satisfactory fit to the data at 
all frequencies, the top panel shows a relatively good fit to 
the higher frequencies and earlier times. This model requires 
a slightly higher electron temperature than typically assumed, 
7’ c j ~ 4 x 10 4 K (although consistent with observations of some 
novae; e.g., Eyres et al. 1996), to fit the high flux densities and 
relatively early peak at 41/45 GHz. It also requires a massive 
ejection (M e j & 1 x 1 0 4 Mq) to achieve the high measured 
radio flux densities. 

Such a large ejecta mass is typical of classical novae, but 
contradicts optical and X-ray observations which imply that the 
nova in V407 Cyg was a relatively low mass ejection, more 
typical of recurrent novae than classical novae. The optical light 
curve declined by three magnitudes in the V band over 24 days 
(Munari et al. 2011). According to nova models, such a rapid 
decline should only occur if M e j < 10~ 6 M q (Yaron et al. 2005). 
In addition, Mikolajewska (2010) points out that V407 Cyg’s 
light curve matches the light curve of RS Oph’s 2006 nova 
event remarkably well; estimates of the ejecta mass for RS Oph 
range across ~10~ 7 to 1 0 6 M 0 (O’Brien etal. 1992; Sokoloski 
et al. 2006; Drake et al. 2009; Eyres et al. 2009; Vaytet et al. 

2011) . Additional support for a low ejecta mass is lent by a 
soft X-ray component, observed in V407 Cyg between Days 
~ 1 5-50 and identified as thermal emission from residual nuclear 
burning on the WD surface (Schwarz et al. 2011; Nelson et al. 

2012) . The speed with which the ejecta become optically thin is 
also confirmed by optical spectroscopy, where high-ionization 
lines (e.g., [Fe xi]) first appear around Day 7 and subsequently 
brighten (Munari et al. 2010). The rapid “turn-on” time of soft 
X-ray emission and coronal lines imply a low ejecta mass of 
~ 1 0 6 to I 0~ 7 M q , assuming an ejecta velocity of 3200 km s“ 1 
(Schwarz et al. 2011). 

A small ejecta mass of M e j = 1 x 1 (L 6 M 0 yields the model 
shown in the bottom panel of Figure 4. This model underpredicts 
the radio flux densities by two orders of magnitude. The model 
light curves also evolve extremely quickly, with the predicted 
5 GHz curve peaking ~20 days after the start of outburst rather 
than at Day ~ 100 as observed. 

Even if the ejection in V407 Cyg was massive, serious 
problems remain in fitting the radio spectral energy distribution 
with the model pictured in the top panel of Figure 4. The 
figure illustrates this point; the model matches the 45 GHz 
light curve rather well up until Day ~ 1 60, but offsets between 
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Figure 4. Models for the multi-frequency radio emission produced by thermal expanding nova ejecta are superimposed as dotted lines on our VLA light curves for 
V407 Cyg. The models assume a “Hubble flow” distribution of ejecta, a maximum velocity of 3200 km s -1 , a distance of 3.0 kpc, and a ratio of Umax/fmin = 0.1. 
The top panel shows a model for a massive ejection that fits the data relatively well (M e j = 1 x 1 0 4 M e , T e j = 4 x 10 4 K), while the bottom panel shows a model 
consistent with multi-wavelength observations of V407 Cyg and expectations for recurrent novae (M e j = 1 x 1 0 6 Mg, 7 C| — 1 x 10 4 K). Note that the vertical scales 
of the top and bottom plots differ. 

(A color version of this figure is available in the online journal.) 


data and models grow toward lower frequencies. At 1.8 GHz, 
there is roughly an order of magnitude discrepancy between the 
observed and model light curves. This discrepancy is caused by 
the model prediction of an optically thick thermal spectrum at 
early times (a = 2), which is much steeper than the a = 0.6-1 .0 
that we observe (Table 2). The model also predicts that, once the 
thermal emission becomes optically thin, it will rapidly fade. 
In contrast, we have observed a late-time plateau in the light 
curve of V407 Cyg with a duration of >300 days (continuing 
through present) during the time when the spectral index implies 
optically thin emission. 

The simple thermal model assumes that the ejecta are not 
clumpy, but smooth in density with a volume filling factor of 
/= 1 . However, observations of novae reveal that the ejecta can 

be quite clumpy, with / = 1 0 5 — 1 0 1 (e.g., Saizar & Feiiand 

1994; Mason et al. 2005; Ederoclite et al. 2006; Shara et al. 
2012). In a simple model, for isothermal ejecta with a single 
filling factor describing all times and radii, dumpiness simply 
increases the radio luminosity by a constant factor over the 
uniform-density case, with no changes to the time evolution or 
radio spectrum (Abbott et al. 1981; Nugis et al. 1998; Hey wood 
2004). This basic case is likely an oversimplification, and more 
complex modeling of the radio emission from multi-temperature 


and time-variable clumped ejecta is needed. However, it is 
difficult to imagine that the addition of such complications 
could slow the evolution of the low-mass model light curve 
by the factor of ~4 needed to match the data. Therefore, while 
clumping might serve to increase the radio luminosity from our 
M e j = 1 x 1 0 () M q model to bring it into better alignment 
with observations, clumping cannot account for the dramatic 
discrepancies in light curve timescales and spectral indices. 

Of course, the nova ejecta in V407 Cyg are decelerating 
with time as they interact with the wind from the Mira giant 
companion (clearly seen in optical line profiles; Munari et al. 
2011; Shore et al. 2011). In the process of decelerating, the 
outer ejecta will be swept up into a relatively thin and dense 
shell; the radio emission will therefore remain optically thick 
longer than in the constant expansion model, and then transition 
more quickly to an optically thin state. Such an evolution 
would only exacerbate the disagreement we already see between 
the data and models of expanding nova ejecta — namely, an 
observed spectral index at early times which is much flatter than 
a — 2 and a more-gradual-than-expected rise and decline in the 
observed light curve. This consideration serves as a reminder 
that significant circumbinary material surrounds the nova, which 
may itself contribute to the radio emission (see the next section). 
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Figure 5. Setup for our model of V407 Cyg. The size of the blue and red circles, denoting the WD and Mira giant, respectively, are larger than to scale for clarity. 
A wind emanates from the Mira giant, depicted as an orange radial gradient. Along a given polar angle 0, the radial density profile of the circumstellar material is 
approximated as independent of azimuthal angle <p. The radial density profiles depend on the separation between the WD and Mira giant (a) and the wind mass-loss 
parameters M w and v w . Two such density profiles are plotted, for 6 = ±20° and 6 = ± 170°, assuming M w = 2 x 10 6 Mq yr 1 , v w = lOkms -1 , and a = 17 AU. 
We calculate r = 1 surfaces and determine their radii from the WD at each of an array of impact parameters (j), assuming that the system is being viewed in the 
orbital plane, with the Mira giant between the viewer and the WD (from the right, in this figure). 

(A color version of this figure is available in the online journal.) 


The required ejecta mass shrinks, but remains quite high 
(~ 4 x 10~ 5 M 0 ) if the distance is small, ~1.5 kpc. Independent 
of distance and clumping, a model with M e j < I O ' 6 M G peaks 
far too quickly and declines too rapidly, compared with our 
observed light curve for V407 Cyg. In view of the poor fits 
to expanding thermal ejecta and the well-established presence 
of circumbinary material around V407 Cyg, we now proceed 
to consider other physical explanations for the radio emission 
from V407 Cyg’s 2010 nova explosion. 

5. RADIO EMISSION FROM AN IONIZED MIRA WIND, 
SHOCKED FROM WITHIN 

We now consider the contribution of the Mira wind to the radio 
emission from V407 Cyg. The ionized envelopes surrounding 
symbiotic stars emit thermal radio emission; their luminosities 
scale with the wind density and the degree of ionization, but span 
the range of radio luminosities observed from V407 Cyg after 
the 2010 nova (Seaquist & Taylor 1990). We show that the Mira 
wind in V407 Cyg can produce the bulk of the radio luminosity, 
as it becomes ionized by the hot WD. The radio luminosity 
grows while the wind is becoming more ionized, and declines 
as this ionized wind gradually recombines and becomes swept 
up by the nova shock. 

5.1. A Quantitative Model 


medium surrounding the WD in V407 Cyg. We solve for the 
ionization of the wind as a function of radius, polar angle from 
the WD, and time. We compare the location of the ionization 
front with the location of the nova shock front (also a function of 
polar angle from the WD and time), to estimate the portion of the 
wind that has been shock-heated and no longer emits in the radio 
band. Then, assuming a viewing direction, we can calculate the 
radio flux density and model the radio light curve. A similar 
mechanism has been suggested by Kenny & Taylor (2005) to 
explain the radio flux densities and images of symbiotic stars 
in outburst; they model these systems as an ionized, fast, hot 
wind emanating from the WD, colliding with a slower, cooler, 
partially neutral wind from the red giant. 

V407 Cyg is an inherently asymmetric system; photons 
emanating from the WD in the direction of the giant encounter 
significantly denser material that those traveling away from the 
red giant. Our model accounts for this asymmetry by splitting the 
environment of the binary into slices of polar angle (6) centered 
on the WD position (Figure 5). Assuming spherical symmetry 
around the Mira giant, a given slice can be described with a 
single radial profile of circumstellar density, independent of 
azimuthal angle (cf>) and time. The density profile in a given polar 
angle slice depends only on the binary separation a, the mass- 
loss rate of the wind M w , and the wind velocity v w . Relative 
to the Mira giant, the density of the spherically symmetric 
circumbinary material is 


Here, we expand on the simple model introduced in Nelson 
et al. (2012), which was originally intended to describe the 
evolution of the nova blast wave in the asymmetric circumbinary 


«CSm(Trg) = 


M w 

4nnmnv w 


,-2 

'RG’ 


(1) 
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where «csm is the number density of particles in the wind at a 
radius r R( -, from the Mira giant, // is the mean molecular weight, 
and »!h is the hydrogen mass. Throughout this paper, we quote 
M w assuming v w — 10 km s _1 , as estimated for V407 Cyg from 
optical PCygni line profiles (Abdo et al. 2010). Relative to 
the WD, the circumbinary medium is asymmetric as a function 
of 9: 


»csm(Avd) = 


M w 

Anixm n v w 


x [(a — r WD cost?) 2 + (r WD sin0) 2 ] , (2) 


where 9 = 0 is the direction pointing from the WD toward 
the Mira giant and rwD is the distance from the WD. Radial 
density profiles for two different values of 9 (one toward the 
Mira giant and the other away) are plotted in Figure 5 . Regions 
of the Mira wind away from the giant will become ionized more 
quickly than toward the giant (illustrated by comparison of our 
model’s density profiles for the 9 = ±170° and 0 = ±20° 
directions, shown in Figure 5). Also, the nova blast wave will 
experience less deceleration in the 9 = ±170° direction than in 
the 9 = ±20° direction. 


5.1.1. Photoionization Model 


We assume that, at the time of nova outburst, the Mira wind 
is completely neutral. We then input a time-variable flux of 
ionizing photons ( N u (t ), in units of photons s -1 ), assumed to 
be emanating from the WD as a result of nuclear burning on 
its surface. Our selections for this important input are described 
more below, in Section 5.1.4. The fraction of the ionizing flux 
that enters a given polar angle slice centered at 0 is 


1 

f\ ol — -COS 



1 

-cos 

2 



(3) 


The width of the polar angle slice is Ad, and set to 1 deg in our 
model (for a total of 180 slices). 

At a given time step, a given polar angle slice is split into radial 
bins of width A r; starting at small twd and traveling outward, 
we calculate how much ionizing flux reaches a given radial bin 
and the ionization fraction x(r, 0, t ) in that bin. The ionizing 
flux at the WD photosphere is S u (r = 0 ,9,t) = f m \N u (t). 
The ionizing flux reaching radial bin (r. 0, t ) is S u (r — 1,0, f), 
while the ionizing flux escaping bin (r, 0, t ), after ionizing the 
neutral atoms in this bin, is S u (r, 6, t). The relationship between 
incident and emerging radiation is 


S u (r, e, t) = S u (r - 1, 6, t ) 

" 47T [1 - x(r, 9, t - 1 )]n csm (r, 0)r 2 lD Ar f m \(9) ' 

At 

- [47 r x(r , 9, t - 1 ) 2 n zsm (r, 9) 2 a (2) % D A r f m \(9)]- (4) 

All parameters are in cgs units. The second term on the right- 
hand side represents the number of neutral atoms in this radial 
bin available for absorbing the ionizing flux. The third term is 
the rate of recombinations in this bin — newly neutral atoms that 
will also absorb S u (r — 1, 9, t). Here, cq 2 ) is the recombination 
coefficient of hydrogen, excluding recombinations to the n = 1 
level. At a temperature of 10 4 K, a. {2) = 2.76 x 10~ 13 cm 3 s -1 
(Spitzer 1978). This value is conservative, as recombination 
plays a relatively weak role in the evolution of radio emission 
from V407 Cyg; hotter assumed temperatures would give lower 
values of cq 2 ) and even weaker effects from recombination. 


We assume that the ionization front is sharp; all atoms in a 
bin must become ionized ( x(r,9,t ) = 1) before the residual 
ionizing flux is allowed to escape to the next radial bin. In 
support of this assumption, we only include photons with 
energies between the Lyman limit and five times the Lyman 
limit (13.6-66.5 eV) in the count of ionizing photons. The 
hardest photons are excluded due to their small cross sections 
for ionizing hydrogen; they presumably serve to ionize metals, 
although we include no treatment of metals in our simple model. 
If the ionizing flux reaching a bin is not sufficient to ionize all of 
the atoms in this bin (in other words, if S u (r, 9, t) in Equation (4) 
is negative), the fractional ionization of this bin is calculated as 


x(r, 9, 1 ) = x(r, 9,t — 1) + 


S u (r — 1 ,9,t)At 


An n csm (r , 9) r^ D A r f m \{9) 
An x{r, 9,t— l) 2 n csm (r, 9) 2 a (2) Ar / vo i(0) At 
An n csm (r , 9) r 2 , D Ar f m \(9) 


.(5) 


This formulation properly accounts for those time steps and bins 
which are not completely neutral ( x > 0) and are not reached 
by any ionizing flux; these will slowly recombine in the absence 
of an ionizing radiation field. 


5. 1.2. Nova Shock Front Model 

Because 10 7 K ionized gas has very different radio-emitting 
properties than 10 4 K gas, we must also track the location of 
the nova shock as a function of time and 9. Increasing the 
electron temperature by three orders of magnitude decreases the 
radio optical depth by four orders of magnitude (Equation (6)). 
The nova blast wave will sweep up circumbinary material 
and heat it from radio-emitting temperatures to X-ray-emitting 
temperatures. 

Our model of the evolving nova blast wave is as described in 
Nelson et al. (2012). We input M w , v w , a, M e j, and initial blast 
wave velocity. The model assumes that the ejecta are expanding 
in a thin shell at a single velocity initially. Each polar angle ring 
then has an independent evolution of radius and velocity that 
depends on the density profile it encounters. It decelerates as it 
sweeps up material, conserving kinetic energy in a Sedov-like 
evolution. 


5.1.3. Calculation of Radio Flux Densities 

The luminosity of thermal radio emission depends critically 
on the emission measure (EM), here defined as f n 2 dl, where 
n e is the electron density of the absorbing material and dl is 
the path length through this material. The optical depth r v at an 
observing frequency v along a given line of sight is 


t„ = 8.235 x 1(T 2 



EM \ 
cm' 6 pc ) 


(6) 


(Seaquist & Bode 2008). The flux density S v produced by this 
line of sight is 



= 3.2 x 10~ 36 






(7) 


Here, A is the projected area of this line of sight, D is the distance 
to the emitting source, and T e is its electron temperature. 
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In our model of V407 Cyg described above, we can calculate 
S v if we assume a viewing angle; optical spectroscopy observa- 
tionally constrains this angle for the 2010 outburst of V407 Cyg. 
Munari et al. (201 1) and Shore et al. (2012) point out that nar- 
row absorption lines of hydrogen and metals are present in the 
spectrum of V407 Cyg following the nova outburst, implying 
the presence of neutral circumstellar material. Although these 
lines fade with time, they never completely disappear, implying 
that the wind never becomes completely ionized. Therefore, the 
nova was probably observed while the giant was in front of the 
WD along our line of sight (Shore et al. 2012). Our best guess 
at the orbital period of V407 Cyg (43 years; Munari et al. 1990) 
is significantly longer than the timescale of our observations, so 
we assume this viewing angle is constant across our radio light 
curve. We show a schematic picture of the binary in Figure 5, 
which we assume is being viewed from the right-hand side of 
the figure. 

While our model is not currently equipped to provide light 
curves for an arbitrary viewing angle, the radio light curves 
presented here should be relatively unaffected by small changes 
in viewing angle, as the radio luminosity is predominately 
dependent on the projected angular area of warm ionized wind, 
which will weakly depend on viewing angle. This insensitivity 
to viewing angle is illustrated for symbiotic stars in Figure 6 
of Seaquist et al. (1984), where, for winds that are significantly 
ionized, the measured flux densities and spectral indices change 
negligibly with a 90° change in viewing angle. 

We assume an electron temperature and construct an array 
of impact parameters (values of z, separated by Az), along the 
axis perpendicular to our line of sight (Figure 5). For a given 
impact parameter, we find all bins of ( r , 9 ) that are intercepted 
by this line of sight, and sum them up. The EM for each bin is 
calculated as 


EM(r, 9, t) = 


x(r, 9 , t) 2 n csm (r, 9) 2 
3.086 x 10 18 

z 


tan (0 + f) tan(<9-f) 


cm 6 pc. 


( 8 ) 


However, a bin is only considered as contributing to the 
radio-emitting EM if it has not been shocked by the nova 
blast wave. If a bin has been swept over by the nova shock, 
calculated as described in Section 5.1.2 and Nelson et al. 
(2012), it will typically be heated to 10 6 — 10 8 K, and will 
therefore be optically thin and faint at radio wavelengths; we set 
EM(r, 0 . t) — 0 in this case. We calculate A by assuming that 
the emitting area at a given impact parameter is ring-shaped: 
A = 7t((z + (Az/2)) 2 - (Z - (Az/2)) 2 ). 

To understand the evolution of the optically thick radio 
emission, we find it instructive to calculate surfaces where the 
optical depth reaches a value of unity, as a function of time. 
Radio observations at a given frequency v will be insensitive 
to conditions within the r v = 1 surface. The size of a r„ = 1 
surface depends on the path length through and density of the 
ionized gas. To reach r v — 1, an EM of ionized material must 
be obtained equivalent to 



Figure 6. Models of V407 Cyg, with the binary located at the origin (small 
red and blue dots; not to scale) and the nova expanding into a Mira wind of 
mass-loss rate M w = 2 x 10 -6 Mq yr _1 , wind velocity v w = lOkrns -1 , and 
wind temperature T e , w = 1 x 10 4 K (orange gradient), as pictured in Figure 5. 
We show the environment of the binary at four different epochs: 1, 25, 200, and 
800 days after outburst. In each epoch, the portion of the wind that is ionized 
is filled with white and outlined in blue. The extent of the nova shock front is 
outlined in green (too small to be visible here on Days 1 and 25). Surfaces where 
optical depth r = 1 is reached are shown as black lines, and are depicted at three 
different frequencies: 1 .4 GHz (dashed), 5.0 GHz (solid), and 25.0 GHz (dotted). 
These surfaces grow as the wind becomes more ionized, but are destroyed when 
the nova shock sweeps over them. 

(A color version of this figure is available in the online journal.) 


given impact parameter, we proceed from the bin closest to us 
along our line of sight, located far out in the Mira wind, inward 
toward the binary, and potentially beyond it, summing up the 
EM in bins until EM T is reached. The bin where EM r is obtained 
marks a location in (r. 9) of the r = I surface. This process is 
iterated for all impact parameters to obtain an array of such 
locations and map the topology of the r = 1 surface (plotted as 
black lines in Figure 6). 

In Figure 7, we plot the evolution of the radius of the r = 1 
surface at three different frequencies, measured perpendicular 
to our line of sight. The radius of the r = 1 surface is defined 
as the maximum impact parameter where r = 1 is reached. 
In some models, the x = 1 surface does not always look like 
a filled circle, but displays a ring-like geometry at later times, 
because EM r is never reached at small values of z (e.g., Day 
800 in Figure 6). Ring-shaped r = 1 surfaces are produced by 
the nova blast eating away at the thermal wind from the inside. 
They are depicted in Figure 7 as having both a minimum and 
maximum radius (e.g., the 1.4 GHz curve in Figure 7, between 
Days 500 and 1000). 

5.1.4. Assumed Physical Parameters 


em t \ 

cm~ 6 pc ) 


12.1 



In exploring the parameter space of V407 Cyg, we hold most 
(9) parameters fixed and focus on varying the wind mass-loss rate, 
M w . The wind velocity is assumed to be 10 km s -1 in our 
discussion of M w . The binary separation is held constant at 
17 AU. In calculating the evolution of the nova blast wave, we 


To calculate the shape and size of the r„ = 1 surface at a given 
time and frequency, we again step through values of z. At a 
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Figure 7. (a) For the same model parameters as depicted in Figure 6, we plot the radius of the r = 1 surfaces, perpendicular to our line of sight, as a function of time 
for 1.4 GHz (gold), 5.0 GHz (red), and 25.0 GHz (blue). The maximum radius of the nova blast wave is also shown as a solid black line, (b) The assumed evolution of 
the WD blackbody temperature, (c) The assumed evolution of the WD blackbody luminosity, (d) The ionizing photon flux produced by the assumed WD temperature 
and luminosity evolution. 

(A color version of this figure is available in the online journal.) 


assume an initial shock velocity of 3200 km s -1 (Munari et al. 
2011; Shore et al. 2011) and an ejecta mass of 10 -6 M e . 

We assume a wind electron temperature of T e w = I x 10 4 K, 
but the calculated radio emission is relatively insensitive to T e w . 
In Figure 8, we plot the same model, one with T e w = 1 x 10 4 K 
(solid lines) and the other with T ew = 3 x 10 4 K (dashed lines). 
At most points in the light curve, the models differ by less than 
a factor of ~2. 

An essential input for our models is the ionizing flux that 
emanates from the nuclear-burning WD, and how this flux varies 
with time. We assume the WD has a blackbody spectrum; theory 
and X-ray observations of V407 Cyg inform our assumptions of 
how the blackbody temperature and luminosity vary with time. 
We plot the assumed evolution of the hot WD in the bottom 
three panels of Figure 7. 

At the onset of the thermonuclear runaway, a brief “flash" 
of ionizing radiation is released, before the ejecta has time 
to expand and reprocess the hot WD spectrum to a puffed- 
up ~ 10,000 K blackbody. We assume this flash lasts 1 hr, has 
a luminosity of 1.2 x 10 39 erg s -1 , and a temperature of 10 5 K 
(Prialnik 1986). Subsequently, the ionizing flux dramatically 


diminishes, until the ejecta have expanded sufficiently to be- 
come optically thin to the hot WD once more. 

There is evidence for significant ionizing flux in the Swift 
X-ray light curve from a soft X-ray component present between 
Days ~ 15— 50 (Shore et al. 2011; Schwarz et al. 2011; Nelson 
et al. 2012). The soft X-ray component is likely attributable 
to the hot WD, which would then be the primary source of 
ionizing flux in the V407 Cyg 2010 nova. Note that the hard 
X-ray component (attributed to forward-shocked material in 
the nova blast wave) is observed to be brighter than the soft 
X-ray component in the Swift light curve, but this is only an 
observational effect produced by the severe absorption of the 
soft X-rays by circumstellar material. Because of the strong 
absorption, it is difficult to precisely determine the luminosity 
of the hot WD from the Swift data (Schwarz et al. 2011; Nelson 
et al. 2012). We therefore assume a WD plateau luminosity of 
1 x 10 38 ergs -1 , as consistently implied by models of nuclear 
burning on WDs during nova events (e.g., Sala & Hernanz 
2005). We do use the Swift measurements to constrain the rising 
temperature of the WD up until a peak temperature on Day 39. 
The rise is fit roughly with 7 WD = 2 x 10 4 1 0 96 K, where t is 
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Time Since Outburst (Days) 


Figure 8. Light curves for V407 Cyg at three frequencies: 1.4 GHz (gold), 
5.0 GHz (red), and 25.0 GHz (blue). The diamond-shaped points connected by 
dotted lines are our VLA data. Models assuming M w = 2 x 1 0 6 Mr. yr 1 and 
T e w = 1 x 10 4 K (as in Figure 6) are shown as solid lines. Dashed lines are 
models for M w = 2 x 10 6 Mg yr -1 and T ew = 3 x 10 4 K. 

(A color version of this figure is available in the online journal.) 

the time since outburst in units of days. Note that this fit is 
extrapolated back to the time of onset of outburst, and implies 
that the ionizing flux from the WD is peaking just a few days 
after outburst (Figure 7). This early time evolution is highly 
uncertain and unconstrained by data, and yet unfortunately it is 
the primary determinant of the early time radio light curve. After 
the temperature peak on Day 39, we assume the luminosity of 
the WD declines from its plateau value as L WD oc r” 1 14 , with a 
roughly constant radius so 7 wd oc f~ 0 - 29 (Prialnik 1986). 

5.2. Interpreting the Light Curve with Our Model 

The ionization state of the Mira wind varies in the aftermath 
of the nova explosion, as evidenced by optical spectroscopy. 
Narrow Ha emission lines are seen at early times (just two 
days after outburst), consistent with an ionized wind (Munari 
et al. 201 1). In the first months after explosion, narrow PCygni 
profiles were observed for low-ionization metal lines (Shore 
et al. 2012). These PCygni metal lines disappeared by Day 70, 
while the Balmer line absorption weakened over this period 
(Shore et al. 2012). Shore et al. attribute this to increasing 
ionization in the Mira wind over the months following outburst. 
This increase in ionization of the wind is not surprising, given 
observations of a hot WD between Days ~ 15-50 (Section 5. 1 .4). 
It is worth noting that some Ha absorption was sustained at all 
times, implying that the wind never becomes completely ionized 
(Munari et al. 2011). 

In our VLA data, the ionized wind of the Mira giant radiates 
as free-free emission, just as detected for other symbiotic stars 
in quiescence (e.g., Seaquist & Taylor 1990). The spectral index 
at early times (a = 0. 6-1.0) is consistent with the theoretical 
expectation for a spherical p oc r~ 2 wind (a = 0.6; Olnon 
1975; Panagia & Felli 1975; Wright & Barlow 1975) and 
spectral index measurements for symbiotic stars (a ranging 
from— 0.1 to 1.4, with a median of 0.8; Seaquist & Taylor 1990). 
However, V407 Cyg is more complex than the typical symbiotic 


star because the radio flux densities are changing rapidly with 
time. Our model shows that the 1.4-25 GHz light curves rise at 
early times because the wind becomes more ionized (Figures 6 
and 8). In other words, the X parameter, as defined in Seaquist 
et al. (1984) to parameterize the fractional ionization of the 
wind, is increasing with time. Seaquist et al. show that as 
X increases, the spectral index of the partially optically thick 
emission flattens — consistent with the decrease in a from 1 .0 to 
0.7 between Days 12 and 39. 

As X increases and the “cone” of the neutral wind decreases in 
extent, the ionized EM at the outer extents of the wind increases. 
The radio flux density thereby brightens and the r v = 1 surfaces 
grow in diameter (e.g., compare Days 1 and 25 in Figure 6). We 
stress that our radio light curve cannot probe conditions internal 
to the x v ss 1 surface. These surfaces are at radii of >10 15 cm, 
significantly larger than the nova blast wave radius (black solid 
line in the top panel of Figure 7). Complex processes may be 
occurring within the r„ = 1 surfaces, like thermal emission from 
expanding nova ejecta (Section 4) or synchrotron emission from 
blast wave interaction (Section 6). However, this emission would 
only become apparent when either the wind becomes optically 
thin or the additional emitting components expand to radii larger 
than the x v ~ 1 surface. 

As can be seen in Figure 8, the model agrees with our 
observed radio light curves as well as can be expected, given 
the significant uncertainties in the ionizing source. The model 
plotted in Figure 8 assumes a Mira wind mass-loss rate of 
M w = 2 x 10 -6 M q yr^ 1 , a distance of 3 kpc, and T ew = 10 4 K. 
It broadly reproduces the observed features of the radio light 
curve. The peak flux densities are roughly matched (to within 
a factor of ~2; compare this with the order-of-magnitude 
discrepancy at lower frequencies in Figure 4), as are the slow 
decline times and the earlier turnover at higher frequency. 

However, the model light curves peak too early, on Day 
~15, relative to the observed maxima of Day ~ 30- 100. Re- 
vising the input ionizing radiation field and its evolution — a 
tweak completely justified given the paucity of observations 
constraining the hot WD at early times — would ease this dis- 
crepancy. We have no constraints on the WD’s luminosity or 
temperature before Day 15 or after Day 50, while its spec- 
trum is peaking in the UV rather than the X-ray — even though 
the WD is likely an important ionizing source at these times. 
We assume that the increase in temperature and luminosity of 
the WD is gradual; extrapolating from the observed temperature 
and luminosity points implies that the ionizing flux is actually 
peaking very early, around Day 4 (bottom panel of Figure 7). 
If, instead, the temperature and luminosity of the WD increased 
suddenly, as observed in RS Oph by Osborne et al. (201 1), the 
predicted rise of the radio light curve would change significantly. 
More detailed modeling of the luminosity and temperature evo- 
lution of the WD is outside the scope of this paper. 

If the mass-loss rate is significantly lower than M w = 
2 x 1 0 6 M ( .j yr _ 1 (solid lines in Figure 9), the flux densities 
never become large enough to match our observations, even if 
the entire wind is ionized. We can understand this by applying 
the standard equation for deriving mass-loss rates of fully ion- 
ized stellar winds from thermal radio flux densities (Wright & 
Barlow 1975): 



12 


The Astrophysical Journal, 761:173 (19pp), 2012 December 20 


Chomiuk et al. 


100.0 


10.0 

js 

5 

a 

Q) 

a 

x 

3 

1.0 


0.1 

10 100 1000 
Time Since Outburst (Days) 

Figure 9. Light curves for V407 Cyg; similar to Figure 8 but showing models 
with different M w . Diamond-shaped points connected by dotted lines are our 
VLA data. Solid lines show a model with M w = 4 x 10 -7 Mq yr -1 . Dashed 
lines are models for M w = 1 x 10 -5 Mq yr -1 ; the 1 .4 GHz model is so faint that 
it falls outside the plot limits. Both sets of model lines assume T etW = 1 x 10 4 K. 
(A color version of this figure is available in the online journal.) 



The Gaunt factor gff is a weak function of the wind electron 
temperature: 


/ 7,3/2 \ 

= 9.77 I —0.17 + 0.130 log^ I . (11) 

We estimate S v from the broad maxima of our 1 .4-25 GHz light 
curves, and find M w ~ (6. 6-7. 7) x 1 0” 7 M G yr -1 for a distance 
of 3.0 kpc and a wind electron temperature of T ew = 10 4 K.This 
estimate is only very weakly dependent on the wind temperature. 
For T ew = 30.000K, M w decreases slightly to (6.0-7. 1) x 
10” 7 M q yr” 1 . These are essentially lower limits on M w . Of 
course, if V407 Cyg is at a smaller distance, the implied M w 
will drop accordingly ( M w ~ 2.5 x 10” 7 M G yr” 1 at 1.5 kpc). 

If the mass-loss rate is significantly higher than M w = 
2 x 10” 6 M Q yr” 1 , as shown by the dashed-line model in 
Figure 9, the hot WD will only be able to ionize a relatively small 
volume. As a result, like in the low M w case, we predict lower- 
than-observed flux densities and small r = 1 surfaces (although, 
in this case the flux density will depend rather sensitively on 
the assumed luminosity and evolution of the ionizing radiation 
field). 

Not only do small r = 1 surfaces produce too little radio 
luminosity, they also lead to radio light curves that decline 
too quickly, because the nova shock front sweeps over them 
at earlier times. For example, in the high M w case plotted in 
Figure 9, the model predicts that the radio emission at 25 GHz 
should be optically thin around the time the observed light curve 
is peaking. 

Radio flux densities eventually decline for two reasons: 
recombination and heating by the nova shock. A good example 
of decline by recombination can be seen in the 1 .4 GHz curve 
in Figure 7; the maximum radius declines from ~ 1 1 x 10 15 cm 
to ~4 x 10 15 because of recombination in the outer parts of the 


wind. At the same time, between Days 600 and 1000, we see 
that the minimum radius of the 1 .4 GHz r = 1 surface increases 
from 0 cm to ~4 x 10 15 cm, until the minimum and maximum 
radii meet and the system becomes completely optically thin to 
1 .4 GHz radiation. 

In our observed light curves, all frequencies stop rising after 
Day ~70 — consistent with the findings of Shore et al. (2012) 
that the wind becomes increasingly ionized up through Day 
~70. Soon thereafter, the observed light curves turn over and 
begin to decline for the higher frequencies (>4 GHz), but the 
1.4/ 1.8 GHz light curves are constant for another ~ 1 00 days 
before declining and becoming optically thin. Correspondingly, 
our model 1.4 GHz curve is quite flat-topped until Day ~300, 
when both recombination and the nova blast begin to signifi- 
cantly diminish the radio EM. 

At late times, the radio flux densities observed at all frequen- 
cies converge to a steady (orperhaps slowly declining) ~3.5 mJy 
(starting around Day 350, continuing through the time of sub- 
mission). The radio spectrum is consistent with optically thin 
thermal emission (a ~ —0. 1 ; Table 2). Such a leveling off of the 
radio light curves is, to date, consistent with our rough model 
for the recombining Mira wind (Figure 8). Occasional monitor- 
ing of V407 Cyg in the upcoming years will show if the radio 
luminosity continues to gradually fade as expected. 

6. WHY NO RADIO SYNCHROTRON EMISSION 
OBSERVED IN V407 CYG? 

From the gamma-ray detection of V407 Cyg, we know 
that particles were accelerated to relativistic speeds by the 
nova shock (Abdo et al. 2010). We might therefore expect 
accompanying radio synchrotron emission. However, at no 
time are our measured spectral indices (Table 2) indicative of 
optically thin synchrotron emission, for which we would expect 
a = —0.5 to —1.0. Here we model the synchrotron emission 
from V407 Cyg to show that bright synchrotron emission was 
likely produced in V407 Cyg, but it was hidden behind the 
absorbing screen of the Mira wind. 

We assume that as the nova shock wave interacts with sur- 
rounding circumstellar material, it converts a constant fraction 
(e) of the post-shock energy density into relativistic electrons, 
such that 


/ N 0 E p dE = e fi /« h «csm vf, (12) 

where all quantities are in cgs units. The blast wave velocity is 
v s and the Mira wind density at the location of the blast wave 
is «csm (as described in Equation (1)). Here, p is the index 
of the electron energy spectrum, assumed to have a power-law 
form. Abdo et al. (2010) estimate p with large uncertainties from 
the gamma-ray spectrum: p — 1 .75//, 4 6 . We can also directly 
measure p if optically thin synchrotron emission is detected, 
as it is in a few other novae. The radio spectral index implies 
p = 2. 1-2.5 in RS Oph (Rupen et al. 2008; Sokoloski et al. 
2008) and p = 2.4 in GK Persei (Seaquist et al. 1989); we 
assume p = 2.3 for V407 Cyg. The minimum energy in the 
integral, 0.5 1 MeV, corresponds to the rest-mass energy of the 
electron and approximates the energy at which electrons become 
relativistic. The efficiency of particle acceleration, e, is assumed 
to be constant with time. Abdo et al. (2010) show that, if the 
gamma-ray emission was solely produced by leptonic processes, 
then e ~ 0.004; if hadronic processes contribute to the gamma- 
ray luminosity, this is an upper limit. 
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Table 3 

Constants for Calculating Synchrotron Emission (Pacholczyk 1970) 


p 

C5 

C6 

2.0 

7.52 x 10" 24 

7.97 X 10" 41 

2.5 

9.68 x 10" 24 

8.10 x 10" 41 

3.0 

1.37 x 10" 23 

8.61 x 10" 41 


The synchrotron luminosity also depends on the magnetic 
field strength, and we consider two physical scenarios for 
its strength and evolution. In the first case, we assume that 
the magnetic field in the forward-shock region is simply the 
magnetic field from the Mira wind, compressed. As suggested 
by Abdo et al. (2010), the magnetic field energy density in the 
wind is expected to be in equipartition with the thermal energy 
density of the wind. We assume that the relevant conditions 
pertain to the wind before it was heated and ionized by the nova 
outburst, so we set the wind temperature to T w o = 700 K as in 
Abdo et al. The magnetic field strength is then simply a function 
of distance from the Mira giant, and is described as 

^comp - = yj32n hcsm ^ T w0- (13) 

Another possibility is that the magnetic field is amplified in 
the post-shock region to be in equipartition with the relativistic 
electrons. This simplification is commonly used to interpret 
the radio emission from supernovae (SNe; e.g., Chevalier & 
Fransson 2006), and was recently utilized by Chugai (2010) for 
outbursts in symbiotic stars. In this case: 


B 


amp — 



M m h «csm 


(14) 


Here, B depends on both the distance from the Mira giant and 
the blast wave velocity. 

The synchrotron flux density of a spherical body of radius R 
at distance D (where both R and D are in units of cm) can be 
found using Equation (1) of Chevalier (1998) to be 


S v = 5.64 x 10“ 22 ^ B~ 1 ' 2 v 5 / 2 


c 6 


x ( 1 — exp 


V — (p+4)/2' 
V\ 


( R/D f mJy. (15) 


Here, C5 and C(, are constants that depend on p as tabulated in 
Pacholczyk (1970), and as listed in Table 3 for convenience. The 
observing frequency is v in units of Hz, and iq is defined as 

/A \2/(P+4) 

V! = 1.25 x 10 19 ijfR c 6 N () j b (p+2)/ (p+ 4) (16) 

We assume that the synchrotron emission is contained in a 
spherical shell between the forward shock and the contact 
discontinuity, typically with a width ~20% of the contact 
discontinuity’s radius (Chevalier 1982). In Equation (16), /is 
the volume filling factor of this emission, so / = 0.42 for a 
width of 20%. 

All the above equations assume spherical symmetry, but we 
know that V407 Cyg is an inherently asymmetric case. Because 
synchrotron luminosity is sensitive to both shock velocity and 
circumstellar density, it is important that we take into account 
this asymmetry in predicting the expected synchrotron signal 
from V407 Cyg. We again use the toy model of Nelson 


et al. (2012) to estimate the evolution of the blast wave in 
this asymmetric medium (see Section 5.1.2). We calculate the 
synchrotron flux density emanating from each polar angle ring 
by using Equation (15) and scaling the flux density by the 
fraction of a spherical shell’s volume that is enclosed in that 
ring (f v oi; Equation (3)). The flux densities from all rings are 
summed together to give the total synchrotron flux at a given 
time, producing the model light curve plotted in Figure 10. 

The light curves in Figure 10 show the evolution of this 
synchrotron model, assuming M e j = 1 0 6 Mq, M w = 2 x 
10 -6 M q yr -1 , v w = lOkms” 1 , € = 0.004, a = 17 AU, an 
initial shock velocity of v s = 3200kms _1 , and a distance of 
3.0 kpc. The left panel assumes a compressed magnetic field 
(Equation (13)), while the right panel assumes an amplified 
magnetic field in equipartition with the relativistic electrons 
(Equation (14)). The synchrotron emission in both models is 
optically thin, with spectral index a = —0.65, so the lower 
frequencies are brighter. The synchrotron flux density decreases 
with time because the blast wave is decelerating and the density 
of the circumstellar material is declining in all directions that 
are significantly contributing to the synchrotron emission. The 
predicted synchrotron emission in the B comp case is much less 
luminous than in the // lnlp case because the magnetic field in the 
forward shock is weaker. 

Both models predict bright synchrotron emission that would 
have been easily detectable with the VLA in the first months 
after outburst, with flux densities ^>1 mJy (a result confirmed 
by the models of Martin & Dubus 2012). However, in reality at 
early times, when the synchrotron luminosity was highest, the 
synchrotron-emitting forward shock suffered free-free absorp- 
tion from the ionized Mira wind (compare the black solid line 
in Figure 7 with the radii of r = 1 surfaces). The large radius of 
the absorbing screen means synchrotron emission is only visible 
at a given frequency once the wind becomes optically thin at 
that frequency. A comparison between synchrotron models and 
observations is most simply carried out on Day ~470, when 
the observed emission from V407 Cyg is consistent with opti- 
cally thin thermal radiation at all observed frequencies. In the 
B comp model, synchrotron emission constitutes a trivial fraction 
of the total flux density on Day 470, even at 1.4 GHz. In the 
/i arnp model, on the other hand, the predicted synchrotron lumi- 
nosities at lower frequencies (<8 GHz) are significantly higher 
than those observed on Day 470, and we might have expected 
to see a clear signature of synchrotron emission in the spectral 
evolution. We stress that no such synchrotron contribution is 
observed in our radio spectra; spectral indices a measured be- 
tween Days 359 and 669 are ~0, as expected for optically thin 
thermal emission ( a = —0.1) and inconsistent with expecta- 
tions for synchrotron emission (a & —0.65). We may be able to 
constrain the magnetic field in the forward shock with our radio 
light curve. 

Of course, our model of the synchrotron emission depends 
quite sensitively on a number of assumed physical parameters. 
Although V407 Cyg is a better-constrained system than many, 
thanks to the community’s wealth of multi-wavelength data, 
there is still room to vary predictions of the synchrotron 
luminosity by orders of magnitude. For example, a decrease 
of the ejecta mass from M e j = 10 -6 M Q to M e] = 10~ 7 M e 
means that the blast wave decelerates much faster, and decreases 
the 1.4 GHz flux density predicted by the B amp model on Day 
470 from 1 1.5 mJy to 0.2 mJy. Because synchrotron luminosity 
depends so sensitively on blast wave velocity, we are likely 
too simplistic in assuming that the ejecta are expelled in a thin 
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Figure 10. Models for synchrotron emission from V407 Cyg, overplotted as dashed lines on our measured multi-frequency light curve (solid lines and diamonds). The 
two panels represent different assumptions about the magnetic field strength in the forward-shocked region: the left panel assumes 5 com p as described in Equation (13), 
while the right panel assumes Z?amp f rom Equation (14). 

(A color version of this figure is available in the online journal.) 


shell with a single velocity; the synchrotron model might be 
improved with more realistic radial distributions of ejecta mass 
and velocity. The B amp model depends on shock velocity more 
sensitively than the B comp model; decreasing M c] to 10” 7 M G 
in the 5 comp case only decrease the 1 .4 GHz flux density from 
0.55 mJy to 0.1 1 mJy. In the fi amp case, the energy densities in 
both the magnetic field and relativistic electrons depend on the 
shock velocity, while in the Z? comp case, the magnetic field is 
independent of shock velocity. 

The synchrotron luminosity also depends sensitively on 
assumptions about the efficiencies of particle acceleration and B 
field amplification in the shock. Although V407 Cyg is unusual 
because e has been directly constrained from the gamma-ray 
observations, this constraint remains a rather uncertain upper 
limit. Additional uncertainties are involved in our assumptions 
about the evolution of e with time and the degree of equipartition 
between relativistic electrons and magnetic field. In the B amp 
case, if e = 0.002 instead of e = 0.004, the 1 .4 GHz flux density 
on Day 470 would have been 3.2 mJy rather than 11.5mJy. 
Assuming the magnetic field and relativistic electrons are in 
equipartition, the magnetic field in the shock on Day 470 must 
have been < 0.009 G to produce less than the observed 1.4 GHz 
flux density. 

While the predictions for synchrotron luminosity remain 
uncertain, reasonable parameters imply that significant syn- 
chrotron emission was produced at early times, but it was hidden 
behind the r = 1 surface of the ionized Mira wind for most of 
the evolution of the nova. Despite rich evidence for particle ac- 
celeration in V407 Cyg, it is not surprising that we did not detect 
synchrotron emission in our VLA data. In contrast, if the ther- 
mal radio emission had been dominated by the nova ejecta, the 
synchrotron emission from the forward shock would not have 


Table 4 

Historical Radio Observations of V407 Cyg 


Date 

(UT) 

Frequency 

(GHz) 

Sv 

(mJy) 

Facility 

Reference 

1982 Jul 13 

4.9 

<0.96 

VLA 

Seaquist et al. (1984) 

1988 Oct 16 

5.0 

<1.72 

BBI a 

Ivison et al. (1991) 

1989 Dec 22 

8.4 

<0.06 

VLA 

Seaquist et al. (1993) 

1992 Oct 10 

270 

<114 

JCMT 

Ivison et al. (1995) 

1993 May 14 

8.4 

1.18 ±0.07 

VLA 

This work 

1996 Aug 18 

1.4 

<2.77 

VLA 

Condon et al. (1998) 


Note. a The Jodrell Bank Lowell-Mark 2 Broad Band Interferometer. 


been absorbed, and would have been detectable — yet another 
reason for dismissing the thermal ejecta scenario investigated in 
Section 4 as the source of radio emission in V407 Cyg. 

7. HISTORICAL RADIO OBSERVATIONS OF V407 Cyg 

Although only radio non-detections of V407 Cyg have been 
reported in the literature for the time period before the 2010 
nova outburst (Table 4), we used unpublished archival data 
from the VLA to determine that the source did in fact produce 
detectable radio emission in 1 993, near the beginning of a period 
of increased activity of the WD (Kolotilov et al. 1998). 

The archival VLA data, obtained on 1993 May 14 as part of 
program AI48 (PI: R. Ivison), were acquired in BnC configura- 
tion at X band with 100 MHz bandwidth; they were calibrated 
using 3C286 and J2007+4029. With 15 minutes on source, the 
1993 data yield a flux density for V 407 Cyg of 1 . 1 8 ± 0 . 07 mJy 
at 8.4 GHz. This detection is in contrast to the non-detection 
obtained four years earlier, when Seaquist et al. (1993) found an 
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upper limit on the flux density of 0.06 mJy at 8.4 GHz. Thus, the 
strength of the radio emission from V407 Cyg varied by more 
than a factor of 20 even before the 2010 nova eruption. 

The radio luminosity of V407 Cyg probably increased around 

1993 because of a flare of ionizing radiation from the WD. 
From spectroscopic monitoring of V407 Cyg over the last two 
decades, it is known that the luminosities of optical emission 
lines vary by factors of ~5 over just a few months. On 1993 May 
29, just two weeks after our radio detection, an optical spectrum 
of V407 Cyg did not show high-excitation signatures of strong 
symbiotic activity (although it did feature Balmer emission, 
which had strengthened since 1991; Kolotilov et al. 1998). By 

1994 August, emission lines of high ionization states (e.g., 
[O hi] 114959, 5007, He n 14686) were visible in the optical 
spectrum for the first time since modern observing began of 
this source in 1984 (Munari et al. 1994; Kolotilov et al. 1998). 
We hypothesize that this active period of V407 Cyg may have 
begun in early 1993 and gradually increased in power. The radio 
appears to precede the strongly symbiotic optical spectrum, but 
note that in spring 1993, the Mira giant was near the maximum 
of its pulsation cycle, and the giant’s light would have diluted 
the optical emission line spectrum. In 1994 August, the giant 
is approaching a pulsational minimum, and the relatively faint 
higher-excitation emission lines would have been more easily 
visible. 

Several years later, in 1998, an outburst was observed in 
V407 Cyg (Shugarov et al. 2007). Unfortunately, no radio light 
curve exists for the interesting period between 1993 and 1998; 
V407 Cyg was not observed again with the VLA (except with 
a shallow NVSS exposure at 1.4 GHz; Condon et al. 1998; 
Table 4) until the 2010 nova. In addition, no information on 
the spectral index exists for the 1993 detection of V407 Cyg, 
so the radio emission mechanism cannot be confirmed. In other 
symbiotic binaries, dramatic increases in radio flux have been 
associated with collimated mass ejection, which can produce 
both thermal and synchrotron radio emission (e.g., CH Cyg, R 
Aqr, and Mira; Taylor et al. 1986; Crocker et al. 2001; Nichols 
et al. 2007; Meaburn et al. 2009). While there is, to date, no 
evidence for fast collimated outflows in optical spectroscopic 
monitoring, V407 Cyg remains a promising candidate site for 
the production of synchrotron emission. The dramatic radio 
variability of V407 Cyg calls for long-term monitoring of 
the radio flux density and spectral index, in cooperation with 
synoptic optical spectroscopy, in order to better understand 
symbiotic activity in this source. 

8. DISCUSSION AND CONCLUSIONS 

Our VLA light curve, spanning frequencies 1 .4-45 GHz and 
covering 17-770 days after outburst, shows a radio evolution 
that is unusually luminous and slowly evolving, compared with 
other novae. Given the postulated low ejecta mass and lack 
of detected synchrotron emission, the most plausible source 
of radio emission from the nova event in V407 Cyg is the 
ionized wind from the Mira giant. The radio emission from 
other symbiotic binaries is explained with a similar mechanism, 
attributing free-free emission to a wind that has been ionized by 
a hot WD (Seaquist et al. 1984). However, the radio emission 
from V407 Cyg is unique due to the relatively rapid rise and 
decline of this radio emission. To achieve a sudden rise in 
the radio light curve requires an abrupt and strong increase 
in the ionizing flux, and the dramatic fading of the radio light 
curve requires a mechanism for disrupting or heating the wind 
on relatively short timescales. We hypothesize that the radio 


light curve declines as the wind slowly recombines on the 
outside and is shock-heated from within by the nova blast. 
The radio light curve provides an estimate of the mass-loss 
rate of Mira, M w « 2 x 10” 6 Mq yr” 1 (for v w — lOkms” 1 
and d = 3 kpc). We place a relatively hard lower limit on the 
mass-loss rate, assuming the wind is completely ionized, of 
M„>7x 10” 7 M q yr -1 . 

8.1. Different Obsen>ational Tracers Imply Discrepant Mira 
Mass-loss Rates in V407 Cyg 

A few difficulties remain for our claim that M w ~ 
I 0” 6 Mq yr -1 in V407 Cyg. In the Suzaku X-ray spectrum of 
V407 Cyg observed on Day 30, non-equilibrium ionization lines 
are observed, implying an ionization age in the forward shock 
of ~2 x 10 10 cm -3 s (Nelson et al. 2012). This relatively small 
ionization age translates to an upper limit on the density of cir- 
cumstellar material at r RG ~ 10 15 cm of «csm 10 5 cm” 3 — or 
a mass-loss rate M w <4x 10” 8 M 0 yr" 1 (an order of mag- 
nitude lower than M w determined from the radio; both assume 
v w = lOkms” 1 ). In addition, studies of the hard X-ray light 
curve can successfully model the data as shocked Mira wind, 
but M w m 1 0” 8 - 1 0" 7 M (J yr” 1 is required (Orlando & Drake 
2012; Nelson et al. 2012). If the wind is too dense, the EM of the 
forward-shocked material will be too high and the ejecta will 
decelerate too rapidly, making the hard X-ray light curve peak 
earlier and brighter than observed (Mukai et al. 2012). 

A similar effect can be seen in the maximum expansion 
velocity as a function of time, measured from the wings of 
He ii 14686 emission line profiles (Shore et al. 2011). If the 
ratio of the wind density to the ejecta mass is too high, the 
ejecta will decelerate faster than observed. We compare the 
redshifted measurements from Shore et al.’s Figure 19 with our 
simple model of the asymmetric blast wave (Nelson et al. 2012). 
The nova material traveling away from the giant companion 
maintains the highest velocity, so we compare this direction with 
the Shore et al. data. We find that the observed velocity evolution 
can be reproduced if, for an ejecta mass of M e j = 1 0” 6 M 0 , the 
wind mass-loss rate is M w ~ 3 x 10 8 M Q yr -1 (again, an order 
of magnitude lower than our radio-determined mass-loss rate). 

On the other hand, the infrared spectral energy distribution of 
V407 Cyg implies M w ~ 8 x 10" 7 M 0 yr” 1 (Kolotilov et al. 
1998), in relatively good agreement with our radio estimate. 

All estimates of M w in V407 Cyg assume that the Mira wind 
is spherically symmetric, which is likely an oversimplification. 
Simulations show that binary interactions lead to a concentration 
of circumbinary material in the orbital plane (Walder et al. 2008; 
Mohamed & Podsiadlowski 2011). Orlando & Drake (2012) 
find that such a density enhancement is critical for matching 
the X-ray light curve of V407 Cyg with a small recurrent-nova- 
like ejecta mass. Radio imaging of V407 Cyg with MERLIN 
and VLA over the weeks to years following the 2010 nova 
show structure in the circumbinary medium (O’Brien 2010; 
T. J. O’Brien et al. 2013, in preparation; A. J. Mioduszewski 
et al. 2013, in preparation). However, the fact that the observed 
radio luminosity is primarily dependent on the area of r = 1 
surfaces implies that we cannot significantly decrease the radio 
estimate of M w by summoning large-scale structure in the wind. 
For example, consider the case of a wind funneling most of its 
mass into a disk-like structure; if the disk is viewed edge-on, 
the r = 1 surface will be relatively small and the radio emission 
will be faint (compared with a spherical wind). The material 
in the disk will be dense, implying that, if we view it face 
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on, it will only be ionized out to relatively small radius (again 
producing small r = 1 surfaces) and will recombine quickly. 
For these reasons, the most effective way to produce large radio 
luminosities is to distribute the wind roughly spherically. 

Small-scale “clumps” in the Mira wind can boost the EM 
of the wind and lead to overestimates of M w from radio data 
(Abbott et al. 1981; Nugis et al. 1998). Small-scale clumping, 
on the other hand, should not affect the IR, X-ray, or Ha 
dynamical estimates of M w . SiO masers are detected in the 
envelope surrounding V407 Cyg (Deguchi et al. 2005, 2011), a 
phenomenon that likely requires an inhomogeneous wind (Gray 
et al. 1998). In addition, the X-ray and radio light curves of 
the 2010 nova show rapid variations (Figure 1; Nelson et al. 
2012), which might be produced by structure in the wind. This 
variety of multi-wavelength diagnostics probe different radii 
in the wind (~10 14 cm from the Mira for masers; ~ 1 0 15 cm 
for the X-ray light curve; ~10 15 -10 16 cm for the radio light 
curve), so we see evidence for inhomogeneities at many radii 
in the wind of V407 Cyg. Therefore, it is possible that our 
radio measurement of M w is an overestimate, and it might be 
brought into better alignment with the X-ray and Ha dynamical 
estimates by accounting for clumping. However, if the filling 
factor of the wind is significantly less than unity, the ionization 
state of the circumbinary material could differ significantly from 
what we predict in Section 5, as it would recombine much more 
quickly. 

Observed mass-loss rates of Mira giants span M w = few x 
10” 8 -few x 1CT 4 M q yr” 1 , with an average ( M w ) = 6 x 
10” 7 M q yr” 1 (Knapp & Morris 1985; although this sample 
is not complete for the lowest M w stars). The Mira giant in 
V407 Cyg has a very long period and is oxygen-rich, so it is more 
typical of OH/IR stars than Miras (Munari et al. 1990). OH/IR 
sources typically show very dense winds with ~10” 3 M e yr -1 
(e.g., Baud & Habing 1983), but V407 Cyg is an odd case 
because dust formation in the wind is suppressed by the hot 
WD, and therefore the wind-driving mechanism may differ from 
typical OH/IR stars (Munari et al. 1990). Still, we might expect 
V407 Cyg to occupy the upper end of the Knapp & Morris range, 
rather than the lower end — more consistent with the radio and 
IR mass-loss rates than those implied by the X-ray data and the 
deceleration of the nova ejecta. 

We therefore postulate that the mass-loss rate in V407 Cyg is 
~ 10” 6 Mp yr” 1 . One possible solution for reconciling the radio- 
derived M w with the lower mass-loss rates implied by X-ray 
non-equilibrium ionization. X-ray emission measure, and shock 
deceleration is a circumbinary medium that is concentrated 
in the orbital plane at small radii (r wl) < 10 15 cm) and 
becomes more spherically distributed at larger radii (r WD > 
10 15 cm). Such a non-spherical wind is in rough keeping with 
the simulations of circumbinary material by Walder et al. (2008), 
which show an overdensity in the equatorial plane of a factor 
of ~1 0-1 00 for r W D < a, and a flare in the outer circumbinary 
disk so that the equatorial overdensity at larger radius is only a 
factor of ~2-3 . However, this scenario would require significant 
fine tuning to come into line with the multi-wavelength data 
on V407 Cyg; the X-ray ionization timescale implies a low 
density at ~10 15 cm while the 25 GHz light curve measures 
M w 10” 6 Mq yr” 1 at ~10 15 cm. A sudden radial transition 
would need to occur in the distribution of circumbinary material, 
from a “disky” concentration in the equatorial plane to a near- 
spherical distribution at a radius of ~10 15 cm to reconcile these 
data. The X-ray and radio data may be brought into better 
alignment if the binary separation is significantly larger than 
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17 AU (in which case, the estimated orbital period of 43 years 
from Munari et al. 1990 would be incorrect). A larger orbital 
separation would decrease the density of circumbinary material 
away from the red giant while still allowing for a relatively 
high M w . Our model for the radio light curve would be largely 
unchanged, but it remains unclear if the X-ray light curve could 
be explained with such a wide binary separation (Orlando & 
Drake 2012 kept a fixed at 15.5 AU in their detailed modeling 
of the X-ray light curve). 

8.2. Implications for Radio Emission from Novae 

The physical process that we claim is responsible for the 
radio light curve of V407 Cyg — namely, the time- variable 
ionization state of the circumbinary environment — has never 
before been recognized as the primary source of radio emission 
in a nova outburst. It provides an interesting test case of 
how dense circumstellar material can affect the evolution of a 
radio transient. Similar outbursts of thermal radiation, produced 
by the flash ionization of circumbinary material, might be 
visible in other novae if significant masses of circumbinary 
material are prevalent, as suggested by Williams & Mason 
(2010). Of course, this emission — particularly if it is optically 
thin — might be superimposed on, or drowned out by, other 
emission components, like the thermal radiation from the 
expanding nova ejecta or synchrotron emission from the nova 
blast. The most similar known nova to V407 Cyg, RS Oph, was 
surrounded by a lower density wind; it did not form an absorbing 
screen following the nova flash, and therefore, synchrotron 
emission was directly detectable from the nova blast (Hjellming 
et al. 1986; Taylor et al. 1989; O’Brien et al. 2006; Rupen et al. 
2008). The recurrent nova T Pyx is surrounded by an Ha+[N ii] 
nebula of mass ~few x 10” 5 M e (Shara et al. 1997; Schaefer 
et al. 2010) which might have been flash-ionized by the 2010 
nova outburst. However, scaling from the Ha+[N n] luminosity 
of this nebula in the months after outburst, we find that the nebula 
produced a negligible < 1 p. Jy of thermal radio emission (Nelson 
et al. 2012). Therefore, V407 Cyg remains the only nova to date 
that shows a radio light curve dominated by the temperature and 
ionization structure of its circumbinary medium. 

8.3. Implications for the Progenitors of SNe la 

The small nova ejecta mass postulated for V407 Cyg implies 
that the WD is massive (>1.25 M Q ; Mikolajewska 2010), 
making this system an intriguing progenitor candidate for an 
SN la. A symbiotic progenitor like V407 Cyg was conclusively 
ruled out in the recent case of SN 201 lfe in M101 (Li et al. 
2011; Horesh et al. 2012; Chomiuk et al. 2012; Margutti et al. 
2012), and has also been ruled out for an additional ~15 SNe 
with deep VLA limits (L. Chomiuk et al. 2013, in preparation). 
However, the sustained elusiveness of SN la progenitors and 
subtle diversity of SN la explosions has led the community to 
consider the possibility that SN la progenitor systems may be 
diverse (e.g., Foley et al. 2012). For example, while most SNe 
la show no evidence for hydrogen in their spectra down to very 
deep limits (Leonard 2007), a very small handful of SNe la 
develop luminous Ha emission lines in their spectra, indicative 
of strong interaction with circumstellar material (Hamuy et al. 
2003; Dilday et al. 2012). In addition, while some SNe la show 
time-variable blueshifted Nai D absorption lines, argued to be 
evidence of circumstellar material, others do not (Patat et al. 
2007, 2011b; Sternberg et al. 2011). 
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Perhaps a fraction of SNe la have symbiotic progenitors, with 
varying delay times between mass loss and SN explosion. For 
those few WDs that explode into a gas-rich medium like present- 
day V407 Cyg, the explosion may resemble the Ha-bright SNe 
la. More commonly, SN la-progenitor WDs are surrounded 
by low-density environments; their immediate environment 
may have been swept clean by nova blasts that swept up 
the giant wind (Patat et al. 2011a; Wood-Vasey & Sokoloski 
2006). However, such a scenario does not describe present-day 
V407 Cyg because, using the classical relation between radio 
flux density and M w (Equation (10)), we derive very similar 
mass-loss rates for V407 Cyg for frequencies spanning the 
range 1.4-25 GHz, which are sensitive to the wind density at 
very different radii (~10 15 -10 16 cm). These observations imply 
that there is no dramatic radial structure in the wind for radii 
>10 15 cm, and the Mira envelope follows a smooth «csm oc r~ 2 
profile. 

Finally, the circumbinary medium around a symbiotic system 
might have time to expand away and merge with the ambient 
interstellar medium, if there is a significant delay between 
accretion and explosion (Justham 201 1 ; Di Stefano et al. 2011). 
V407 Cyg could constitute such a progenitor in ~ 10 5 years time. 
In this case, any blueshifted Nai D absorption lines present in 
the resulting SN la spectrum would be the product of the Mira 
envelope slowly drifting away from its parent binary. 
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